Resonant effects on the two-photon emission from hydrogenic ions 
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A theoretical study the all two-photon transitions from initial bound states with rii = 2, 3 in 
hydrogenic ions is presented. High-precision values of relativistic decay rates for ions with nuclear 
charge in the range 1 < ^ < 92 are obtained through the use of finite basis sets for the Dirac 
equation constructed from B-splines. We also report the spectral (energy) distributions of several 
resonant transitions, which exhibit interesting structures, such as zeroes in the emission spectrum, 
indicating that two-photon emission is strongly suppressed at certain frequencies. We compare two 
different approaches (the Line Profile Approach (LPA) and the QED approach based on the analysis 
of the relativistic two-loop self energy (TLA)) to regularize the resonant contribution to the decay 
rate. Predictions for the pure two-photon contributions obtained in these approaches are found to 
be in a good numerical agreement. 
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I. INTRODUCTION 

Two-photon transition in hydrogen and hydrogenlike 
ions are under investigation since Goeppert-Mayer pre- 
sented her theoretical formalism in 1931 [l|. The early 
interest on these transitions from metastable states of 
hydrogen came mainly from astrophysics which 
was recently revived by Chluba and Sunyaev [J] . Among 
many applications of recent two-photon studies, one can 
cite the determination of the Rydberg constant d, d, 0| j 
measurement of the Lamb-shift 0, testing Bell's in- 
equality 0, as well as various ap plic ations in molecular 
spectroscopy llOll . tissue imaging |ll| and protein struc- 
ture analysis |12l | . Another interest in two-photon transi- 
tions is connected to the stud y o f parity- violation effects 
in H-like and He- like ions (l3l.TT^. The two-photon spec- 
tral distribution has recently been used for precise effi- 
ciency calibration of solid-state X-rays detector as it has 
a known shape for a large distribution of energies [l^ . 

Similar to single-photon processes, two-photon emis- 
sion can be spontaneous or stimulated, whereas two- 
photon absorption is only stimulated. However, since 
each photon carries one unit of angular momentum in the 
dipole approximation, certain transitions between atomic 
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energy levels, forbidden as single-photon processes, are 
allowed as two-photon processes. Another important dis- 
tinction lies in the fact that the emission spectrum of 
spontaneous two-photon transitions is continuous unlike 
the spectrum in a single-photon process. A continuous 
spectrum is possible because energy conservation requires 
only that the sum of both photon energies equals the en- 
ergy of the transition. For the transition 

{ni.,ii) ^ {nf,jf) + hioi + fujJ2., (1) 

where (ni,ji) and {nf,jf) denote the principal quantum 
numbers and total angular momenta of the initial and 
final hydrogenic states, respectively, and hioi and hu!2 
are the energies of each photon, the conservation of the 
energy leads to the condition 

Ef - Ei = huji + huj2, (2) 

where Ei and Ef are the energies of the initial and final 
ionic states, respectively. 

Because of its importance, the 2si/2 lsi/2 two- 
photon transition rate in hydrogen has been calculated 
and discussed many times using different approaches. An 
historical overview from both theoretical and experimen- 
tal point of view can be found in the 1998 article by 
Santos et al 16]. 

Recently, Surzhykov et al ^ITj performed a relativistic 
calculation to study the angular correlations in the two- 
photon decay of hydrogenlike ions and Labzowsky et al 
[ISt evaluated the 2E1 contribution for the 2si/2 lsi/2 
transition and the ElMl and E1E2 contributions for 
the 2pi/2 — > lsi/2 transition using an expression similar 
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to the one obtained by Goldman and Drake [19| in the 
Quantum Electrodynamics ( QED ) framework. Also in 
this framework, Nganso et al [20| carried out the treat- 
ment of the iS* matrix for bound-bound transitions. 

In this work, which uses techniques of a previous one 
[l6l |. we study the two-photon decay of several excited 
states using two approaches to deal with resonances; the 
Line Profile Approach (LPA) [2l[ and the QED approach 
based on the analysis of the relativistic two-loop self en- 
ergy (TLA) to re gula rize the resonant contribution to 
the decay rate 13, 2^- We present calculated values for 
two-photon decay rates obtained with both approaches 
for one-electron ions with a nuclear charge up to 92. This 
article is organized as follows: in Sec. [ill we give a brief 
review of the background theory involved in two-photon 
emission, in Sec. IIIII we present the results obtained in 
this work and the conclusions are presented in Sec. IIVI 



II. 



THEORY OF RELATIVISTIC RADIATIVE 
TRANSITIONS 

A. Two-photon spontaneous emission 

1. General formalism 



The relativistic theory of two-photon transitions is 
given in detail in references [H, [l^ [13l ■ 

In the present work, therefore, we report only the most 
important equations and notations used. 

The basic expression for the differential (in energy of 
one of the photons) rate is, in atomic units. 



dw 
duji 



LU1UJ2 



(2^) 



E 



{f\A*\u) {,.\Al\i) 



{f\Al\u) {u\A*,\^) 

Ey ~ Ei + UJ2 



dVtidVl2, (3) 



where Uj is the frequency and c?f2j is the element of solid 
angle of the j*'' photon, and c is the speed of light. The 
frequencies of the photons are constrained by the energy 
conservation Eq. ([2]). 

For photon plane- wave with propagation vector kj and 
polarization vector Sj {ej -kj—Q), the operators A* in ([3|) 
are given by 



Gk, 



-iki 



^ - Ge- 



iki 



(4) 



where a arc Dirac matrices and G is an arbitrary gauge 
parameter. Among the large variety of possible gauges. 
Grant [1^ showed that there are two values of G which 
are of particular utility because they lead to well-known 
nonrelativistic operators. If G = 0, one has the so 
called Coulomb gauge, or velocity gauge, which leads to 
the dipole velocity form in the nonrelativistic limit. If 
G = [(i + 1) for example, G = \/2 for El tran- 

sitions (L = 1), one obtains a nonrelativistic expression 



which reduces to the dipole length form of the transition 
operator. The two-photon transitioi is g auge invariance 
was studied by Goldman and Drake From the gen- 
eral requirement of gauge invariance the final results must 
be independent of G. 

The index stands for all solutions included the dis- 
crete and both negative and positive energy solutions of 
the Dirac equation. In Eq. ([3]), moreover, \i) — \niKimi), 
W) = \n,yKiym,j) and |/) = \nfKfmf) are the well-known 
solutions of the Dirac Hamiltonian for a single electron, 
where n and m stand for the principal quantum num- 
ber and the one-electron angular momentum projection, 
respectively. The Dirac quantum number k is defined by 



I if j = ^-l/2 

-(^+1) if J =^+1/2 



(5) 



where £ and j are the electron orbital and total angular 
momenta, respectively. 

If the energy of an intermediate state E^, is equal to 
the energy (£"— a;i_2) in the denominators of Eq. ([3]), the 
differential emission rate has a pole or a resonant behav- 
ior at E,y. Physically, this occurs when an intermediate 
virtual state, between the initial and final states, coin- 
cides with a real state so that the two-photon transition 
coincides with the cascade de-excitation process. 

For example, in the 2E1 3si/2 — *■ lsi/2 transition, 
the shape of the frequency distribution presents nar- 
row resonances at energies corresponding to the 3si/2 — > 
2pi/2,3/2 ^ lsi/2 cascade. This effect has been confirmed 
both experimentally [2^ and theoretically [27| . 

The divergent behavior of the resonant denominator 
in Eq. ^ is related to the Green function used in that 
expression, which does not take into account the interac- 
tion between the electron and the vacuum fluctuations of 
the electromagnetic field. The LPA allows to derive the 
following expression for the differential emission [2l|, |2^ , 
that takes partially into account this contribution. 



dw 



LPA 



UJiU!2 



E 



{f\A*2\i^) (H^tN) 



dui (27r)3c2 



K - K + UJ2 



2 

dnidn2, (6) 



where 



with 



1 if is a resonant intermediate state 
" otherwise ' ^ ' 

and {f \ J2e \ ^) I He I ^) ^'^^ electron mean value 

of the self-energy and vacuum polarization operators in 
lowest order for the state respectively. Both the mean 
value of the self-energy and vacuum polarization oper- 
ators have a real part, AE'^, that is a correction to the 
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energy E^. On the other hand, only the self-energy op- 
erator has an imaginary part, ry/2, which is the width 
of the state i/. 

The average decay rate, i.e., the decay rate summed 
over the final ruf and averaged over initial rrii ion mag- 
netic sublevels, can be obtained from Eq. ([6|) as 



dW 



LPA 



duJi 



E 

Z/1 , Ai ,Z/2 , A2 



dW 



LPA 



duJi 



(9) 



where the partial decay rates describing the two-photon 
transitions of a given type (A) and multipolarity (L) are 
given by 



dW 



■LPA 

-Li , Ai ,-^2 , A2 

duJi 



UJ1UJ2 



(2^)3c2(2j- + 1) 



E 



^'"(2,1) + ^'"(1,2) 



; {rc [S'^'^(2,l) 



Im 



^■''^(2,1) 



Re 
Im 



S'^''(l,2) 



■(1,2)]} 



Here we define 



which represents the angular coupling, and 



(10) 
(11) 



^'(2,1) = E 



M 



An 



<A2,-t2) . s -jrj(Al,Ll) , , 
(^2) ^ [UJl) 

r rJi/2 ^»(l)^/(^) 



^' 1 0-1 



with 



4it) 



1 if 4 + ^ + + At = odd 
if + ^ + Lt + At = even 



(12) 



(13) 



£'"'(1,2) is analogously defined. The notation [j, A:,..] 
means (2j + l){2k + 1) . . ., (• • • ) are the 3j symbols and 
{• • • } the 6j symbols. 

The radial matrix elements M^^-^^ in Eq. (fT2| are 
defined by 



L + 1 



1/2 



[(«/-«.)^^+l + (i + l)^L+l] 



L + 1 



1/2 



(14) 



Mff 



2L + 1 



and 



M 



[L(L + 1)] 



G[(2i: + 1) 

-f (k/ - K,;) (/ 

-LI 



(15) 



'L-1 



L+1 ' -"L- 

(i + i)^Z+i] 



-1) 



(16) 



L is the photon angular momentum and A stands for 
the electric (A=l), magnetic (A=0) and the longitudinal 
(A — —1) terms. We used the notation given by Rosner 
and Bhalla j22| for the integrals the if (uj) and (lu). 
The parity selection rules (fT3|) follow from the calculation 
of the reduced matrix elements expressed in Eq. 

We emphasize that the term 7rf(l)7r^(2) in (|12p is not 
given explicitly in the Goldman and Drake article p^ . 
which could lead to some ambiguity in the choice of the 

intermediate states for the evaluation of the S'"'(2, 1) and 

S'"'(l, 2) terms in a generic transition. 

Usually, it is convenient to express the results in terms 
of the electric (E) and magnetic (M) multipole contribu- 
tions. The total decay rate (integrated over the photon 
energy) for a transition in which one photon QiLi and 
one photon 62^2 are emitted, where Qi = E, M stands 
for the electric and magnetic multipole type, respectively, 
is given by 



-LPA 

eiLie2L2 



E ^LuXe,,L 

tLPA 



2,Ae2 



with 



E 



Ae. = -1,1 
Ae, =0 



uJt dW 



L\ ,A0^ ,L2 Ae2 
dhJi 



dtui, 



if 
if 



e^ = E 
e,, = M 



(17) 



(18) 



and ujt is the energy of the two-photon transition, which 
is given, in a.u., by 



(19) 



using Eq. 

Finally, the total spontaneous emission probability per 
unit time for a two-photon transition is obtained by sum- 
ming over all allowed multipole components. 



^ *eiLi,e2L2W^eiLie2-L2, (20) 
all eiLi,02L2 



where 



1 if OiLi 7^62^2 
ieiLi,e2L2 - S 1/2 if OiLi = 62^2 ■ ^ ' 
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The factor 1/2 is included to avoid counting twice each 
pair, when both photons have the same characteristics. 

Another method for dealing with resonances was de- 
veloped by Jentschura and co-workers [H, using a 
procedure based on two-loop self-energy (TLA). They 
obtained an expression similar to Eq. Q for evaluating 
a nonresonant component of the two-photon decay rate, 
given by 



„TLA 



lim Re 



duji 



LO1UJ2 



(2^) 



3^2 



(22) 



The function Sij is given, as in Ref. [23], by 



E 



E,j — Ei 



{f\Al\v) {u\Al\z) 
— Ei + UJ2 — ie 



(23) 



Using this approach one obtains finite results since the 
integration over the frequency uji is displaced by a in- 
finitesimal quantity, e, from the resonance poles, provided 
the limit is not permuted with the integration. 

If one considers a nonresonant transition such as 
2'Si/2 ~* 151/25 then the limit can be permuted with the 
integration and Eq. ([22|l reduces to Eq. ([3|), and both 
approaches gives the same result. 



does not contain a resonant behavior. It is defined as 



,LPA 



(-1) = E 



(25) 



The coefficients and a\ are derived from the Taylor 
expansion of (o^i) around wr: 



Oo = /■^ (^r) = 9' (^r) 




duji 

d 

dioi 



P (^1) 



9^ (^1) 



(26) 



The expressions of the derivatives of the matrix elements 
used to evaluate a\ are presented in the Appendix. 

LPA 

To obtain the decay rate W L2 \2 must add to 
the integral of the smooth function h^^^ the two terms 
()o^^ and evaluated analytically, i.e.. 



W 



-LPA 

-Li , Ai ,-^^2 5 A2 



LPA 



LPA I hLPA 



(27) 



where 



2. Integration method for resonant intermediate states 

For resonant transitions, Eq. ([TO)) produces sharp 
peaks near the resonant frequencies, which requires spe- 
cial attention in the integration over the photon energy 
(jji in Eq. ([TT]) to avoid meaningless results for the total 
decay rate. Near a resonant frequency w;^" , Eq. PH]) can 
be written as 



dW 



LPA 

-Li,Ai,L2,A2 



dull 



p (^1) 



(24) 



where (uji) is a smooth function; the resonant behavior 
is given by the denominator. Consequently, the function 
(wi) can be expanded in a Taylor series around the 
resonant frequency co^ . Notice that the shape in the 
right-hand side of Eq. ((24)) is not a Lorentz profile since 
f^{ijji) depends on uji and so, the peak profile is asym- 
metric. Subtracting the first two terms of the expansion 
on (wi) we obtain a smooth function, h}"^^ which 



LPA 



LPA 



,LPA 



LPA 



duJi 



(28) 



E< 



^1 - - 



jdioi 



E 2Qo 

J ^R 



arctan 



^ R 



,(29) 



■Jo 



TrdlOl 

^1 - ^R - 



r^r/2 



.(30) 



We note that Oq is given approximately (unless we con- 
sider the limit F ^ 0, in which case it is given exactly) 
by 



2tt 



(31) 
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where the term wf^^'' (w) is the decay rate from a initial 
to a final state through the emission of one photon, which 
is given, in a.u., by [25| 



_ f is L j 



c[L] 







M 



Using this result we can write the term 



y- 

3v 



(32) 



(33) 



and identify f}Q^'* as a cascade transition rate contribu- 
tions. 

Applying a similar approach to w^^^ given by Eq. ([^ 
we obtain a smooth function K^^^ as in the LPA. One dif- 
ference between the two approaches is in the term of order 
- (rR/2)', which appears in the denominator of Eq. (j25|) 
and results from considering the infinitesimal quantity e 
finite, i.e., taking the role of a level width (e ^ F). In 
the present evaluation we obtain the function K^^^ by 
replacing Fj^ — > gF^, where g is a parameter that can 
be made arbitrarily small. We thus obtain convergence 
since the difference in li^^^ using g = 1 or g = 10^^ is 
in the fifth digit. For q = 10~^ and q — 10^'^ the differ- 
ence in K^^^ is in the ninth digit. So we conclude that 
using h}^^^ defined in Eq. (|^5|) with q = 10~^, is a good 
approximation for the function h^^^ . Another difference 
between TLA and LPA is the inclusion of radiative cor- 
rections Re [SE] and VP, which for values of Z as high 
as 92 changes the value of h from one approach to an- 
other in the second digit. The major difference between 
the two approaches is in the integral [}o, which in TLA is 
given by 



1 



J. 



(34) 



which comes from the different ways the pole regulariza- 
tion is done. 

Notice that the terms and \]^^^ are related by 



LPA 




TLA 



^ R, 



(35) 



which shows that the difference between and f) J'^^ is 
mainly due to the second term on the left side of Eq. ([551) 
(since F-^ ^ 1) or by the product of one-photon transi- 
tions (cascade process). 

On the other hand, the integral f)i is given in the TLA 
approach by 



TLA 



^<ln 



(36) 



This expression can be obtained from Eq. ([50)1 by taking 



Considering that 



:,TLA 



lTLA 



-TLA 



duJi 



(37) 



the decay rate in the TLA, Wl/^i La given by an 

expression similar to Eq. (|27p . in which the LPA contri- 
butions are replaced by the correspondent TLA contri- 
butions. 

As we will see in Section lllll these differences on the 
sum of [) and [)i do not carry any sizable difference be- 
tween the LPA and TLA methods for low-Z ions, but 
lead to slight discrepancy for heavier systems. 



B. Solution of the Dirac-Fock equation on a 
B-splines basis set 

To make the numerical evaluation we consider that the 
atom, or ion, is enclosed in a finite cavity with a radius 
large enough to get a good approximation of the wave- 
functions, with some suitable set of boundary conditions, 
which allows for discretization of the continua. 

Let us denote by {(/)Jj(r), i = 1, ..,2A^} a set of solu- 
tions of the Dirac-Fock equation, where n is the level, i 
the position of the solution in the set and N the number 
of functions in the basis set. 

For each n value the set is complete, and (/'Jj(r) obeys 
the equation [s^l 



V(r) 



\ dr 



_d_ 

dr 



-2c- 



R 

y(r) 



<Pn{r) - -fnir), (38) 
c 



where the energy i?^ was replaced by el^ — E^^ — mc^ . 
The potential V{r) is given by a Coulomb potential as- 
suming a uniform nuclear charge distribution for a finite 
nucleus and k, is given by Eq. A complete set spans 
both positive and negative solutions. Solutions labeled 
by i = 1, iV describe the continuum < —2mc^ and 
solutions labeled by z = N + 1,...,27V describe bound 
states (the few first ones) and the continuum > 0. 
For practical reasons, such as easy numerical implemen- 
tation, this set of solutions is itself expressed as linear 
combination of another basis set. We have chosen the 
B-splines basis set and we used the derivation of the so- 
lution of equation psp in terms of B-splines described by 
Johnson, Blundell and Sapirstein [3l| . 



III. NUMERICAL RESULTS AND DISCUSSION 

By taking t]^ = Q in Eq. ([7]) we may calculate the 
two-photon decay rates without accounting for radiative 
corrections. In this case, we have verified with respect 
to variations of the gauge parameter (G = for velocity 
gauge and G = y/2 for length gauge), the radius of the 
cavity (i?), and the basis set parameters (the number 
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(ns) and the degree (fc) of the B-spHnes), the stabiHty 
and accuracy to six digits on the calculation of Eq. (fTU]) 
for a nonresonant states and for a frequency loi. 

The parameters used in the calculation of the results 
presented in this work are fc = 9, ns = 60 and i? = 60 
a.u.. The integration over the photon frequency has been 
performed using a 15 points Gauss-Legendre algorithm 
for the nonresonant transitions. 



A. Nonresonant transitions 

For the nonresonant 2si/2 ^ lsi/2 Sind 2pi/2 l'Si/2 
transitions, we use Eq. ([9]) for both the decay rate values 
for various multipole combinations {OiLi, Q2L2) and the 
frequency distribution. 

The most significant multipole combinations included 
in the calculation of the two-photon decay rates of the 
2pi/2 ~^ Is 1/2 transition are presented in Table [H The 
magnitude of the multipole combinations not listed in 
this table are, at least, four orders of magnitude smaller 
than the most significant 

The values given by Labzowsky et al. in were 
obtained using expressions similar to the ones used by 
Goldman and Drake [l^. In Ref. [s^, the results were 
obtained using nonrelativistic Coulomb Green's function, 
which, for high values of Z such as 92 leads to inaccurate 
values. The relative difference between our results and 
the results in Ref. is in the range of 0.1 - 0.4%. We 
observe that, for the three studied Z values, more than 
99% of the total decay rate is due to the multipole con- 
tributions £:iAfl (about 60%) and EIE2 (about 40%). 
The fact that the two multipole combinations ElMl and 
E1E2 give almost the same contribution is somewhat ex- 
pected since Ml and E2 have the same order of magni- 
tude in the decomposition of the photon field [23]. On 
the other hand, a comparison between the listed most 
significant (ElMl) and less {E2E3) significant contri- 
butions reveals that the relative importance of the latter 
increases with Z, being 12, 5 and 4 orders of magnitude 
smaller than the former for Z = 1, Z = 40, and Z = 92, 
respectively. 

In Table nil we report the two-photon total decay rates 
for 2pi/2 ~* lsi/2 and 2si/2 lsi/2 transitions. Enough 
multipoles have been included in the calculation of the 
total 2-photon decay rates to reach an accuracy of six 
digits. The values for the 2si/2 ~* 151/2 transition differ 
slightly from the ones in our previous work (l6j due to 
use of the most recent values of physical constants [s^ , 
such as the fine-structure constant. 

It sould be mentioned that the interest in the transi- 
tion 2pi/2 — > lsi/2i and other two photon forbided tran- 
sitions, is only academic since the transition is suppressed 
by selection rules and this channel is in direct competi- 
tion with an allowed one-photon transition. 

To present the spectral (or frequency) distribution for 
a specific value of Z is convenient to express the results 



in 2lj{y, Z) as suggested by Spitzer and Greenstein 



dW 
dy 



no 



(Za)" ip{y,Z), 



(39) 



where y = oj/ujfi is the fraction of the photon energy 
carried by one of the photons and ujfi is the energy of 
the transition. In case of an even^even (or odd— >odd) 
transition, the major multipole contribution 2E1 scales 
as Z^ and, consequently, n = 6. For a even^ odd (or 
odd^cvcn) transition, both ElMl and E1E2 scale as 
Z». 

In Fig. [Jl the frequency distribution of the multipole 
contributions ElMl, MlEl, E1E2 and £'2£'l for the 
transition 2px/2 ~* ^^1/2 are presented. Although each 
one of these four most significant contributions is asym- 
metric, the sum of each pair {ElMl, MlEl) and {E1E2, 
E2E1) is symmetric around y — 0.5. Therefore, the to- 
tal frequency distribution is also symmetric around the 
y — 0.5 value, as can be seen in Fig. [51 in which we 
also notice the Z dependence of the shape predicted by 
Goldman and Drake [l^ for the 2si/2 



1*1/2 transition. 



Resonant transitions 



After this brief discussion of the nonresonant 2si/2 
lsi/2 and 2pi/2 lsi/2 two-photon transitions, we now 
turn to the evaluation of the differential total decay rates 
for the higher excited ionic states. In Fig. [31 for example, 
we display the spectral distribution for the 3si/2 lsi/2 
transition. We notice several features that are not found 
in the corresponding plot for the 2si/2 ^ lsi/2 tran- 



sition. In particular, the V'3si/2^isi/2 {2/j ^) function 
exhibits sharp peaks, which are due to the 3si/2 
2pi/2,3/2 lsi/2 cascade. Furthermore, we observe that 
at Z = 92 each of the two resonances splits in two due to 
the spin-orbit interaction and the frequency gap in each 
pair is exactly equal to the difference between the states 
2p3/2 and 2pi/2, respectively. In addition, besides the 
zeroes at the endpoints, there are two more minima at 
y = 0.219733 and 0.780267. Such minima were observed 
in two- photon spectra by Tung et al [13, HH] , and they 
were referred to as "transparencies" . In Table IIIII we list 
the transparencies for several two-photon transitions ob- 
tained in this work by other authors. Their relative dif- 
ferences are smaller than 0.01% for Z = 1. To the best 
of our knowledge, there are no published data for other 
Z values. In Fig. [H we plot the transparency frequency, 
ytransp^ of the transition 3si/2 — > 1*1/2 as function of Z. 
We notice that the transparency values scale with Z^ as 
the transition energy. 

In contrast to the 3si/2 1*1/2, the spectral distri- 
bution for the 3^3/2 ~^ 1*1/2 transition, plotted in Fig. 
[5l exhibits only the resonant behavior as mentioned in 
Ref. [34|, which is due to the fine-structure splitting be- 
tween 2pi/2 and 2^3/2 and 3pi/2 and 3c?3/2 states. 

In Fig. [SI we plot the frequency distribution of the 
multipole ElMl contribution for the 2p3/2 1*1/2 tran- 
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sition. Along with the resonances, the shape of the curve 
is similar to the one in Fig. [1] for the l-p^j^ lsi/2 
transition. In the ElMl case, the resonance in the low- 
frequency side occurs when the energy of one of the pho- 
tons is equal to the energy difference -E2P3/2 ~£'2si/2 ; while 
the resonance in the high-frequency side occurs when the 
energy of one of the photons is equal to the energy dif- 
ference £;2p3/2 - ^2pi/2 ■ 

The list of the radiative corrections contributions for 
some states, which were included in Eq. to achieve 
an accuracy of at least six digits are listed in Table IIVI 
The values for the real part of self-energy and vacuum po- 
larization were obtained from the MCDF code developed 
by Desclaux, Indehcato and collaborators [H, 113, Is^ . 
The level width, F^, is equal to the sum of the one- 
photon partial level widths, given by Eq. ([5^ . 

As seen from Eq. p7p . by performing the integration 
of the differential transition probabilities over energy of 
the emitted photon we may finally obtain the total two- 
photon decay rates. Eq. (p7)) shows that these rates can 
be traced back to t) functions. In Table |Vl we list the 
sum of the terms i}^^^ and f)i^^ given by Eq. (^5)1 and 
Eq. ([50]) , required for the evaluation of the decay rates 
for transitions from bound states with = 3 in the 
LPA, including the most relevant multipoles, radiative 
corrections and using q — I. The correspondent values 
obtained in TLA are listed in Table IVII By comparing 
the values in these two tables we conclude that they differ 
less than 0.001% for Z = 1, 2.3% for Z 40 and 10% 
for Z — 92, which shows the importance of the radiative 
effects. 

In Tables I VIII and IVIIIl we list the most relevant 
multipole combinations included in the calculation of 
the two-photon decay rate for the 2p^/2 ~^ ^^1/2 s-nd 
3si/2 ~^ 2sx/2 transitions. 

We notice that for Z = 1 the decay rate values of 
some multipole contributions, such as the ElMl and the 
E1E2, listed in Table IVlTl are similar to the correspon- 
dent ones for the transition 2pi/2 lsi/2- Nevertheless, 
this is not the case for Z = 40 and Z = 92. This is 
due to the fact that the energy separation between 2^3/2 
and 2si/2 increases with Z and, consequently, the decay 
rate contribution from the cascade process also increases. 
This aspect is also evident in Fig. [3 where the multipole 
combination ElMl decay rate Weimi, obtained in the 
LPA and TLA, is plotted as a function of the atomic 
number Z for the 2px/2,3/2 ~^ lsi/2 transitions. 

The resonant behavior of the 2^3/2 — > lsi/2 transi- 
tion is strongly suppressed for low Z values. We no- 
tice that for lower Z values both solid (LPA) and dot 
(TLA) lines have similar values, which is a consequence 
of the fact that nonresonant contribution (related to in- 
tegral of "background"), in both transitions (MlEl in 
Figs. [T] and O is much higher than the cascade term 
(dash line). For higher values of Z, we notice that the 
solid line follows the dash line. This could be explained 
by the different Z scaling of the two contributions. The 
"background" scales as Z^ and the cascade term, given 



by 2p3/2 — > 2si/2 — > lsi/2i scales as Z^'^. The dash-dot 
(decay rate of 2pi^2 ~ 1*1/2) dot lines are almost co- 
incident in low Z region and diverge from about Z = 40, 
which is an evidence of the relativistic effects in the npi/2 
and np3/2- 

In Table HXl we report two-photon total decay rates for 
transitions from initial level with rii = 3, obtained in the 
LPA considering the most relevant multipole combina- 
tions in Eq. (PO)) . The results of Tung et al, presented in 
this table, were calculated using the analytical formulas 
described in Ref. [l^l, which were obtained through the 
so called implicit technique that describes the intermedi- 
ate states by a differential equation. 

We restrict ourselves to list the two-photon decay rates 
obtained in the LPA because in some cases they are very 
different from the TLA ones, when the cascade term in 
Eq. (|55)l dominates. 

One important aspect concerning total decay rates of 
resonant transitions is the calculation of the nonresonant 
decay rate without interference from resonant intermedi- 
ate states. Cresser et al [39^, using a fourth-order pertur- 
bation term development, obtained an expression similar 
to Eq. (35.21) in Ref. [2J| where the sum over the in- 
termediate states considers only the states above the ini- 
tial one, avoiding in this way the resonant denominators, 
and found the value 8.2197 s~^ for the 3si/2 lsi/2 
transition rate. Florescu [i^l, using the same procedure, 
obtained the value 8.22581 s~^. The nonrelativistic limit 
of Eq. ^ in the Coulomb gauge gives the same ex- 
pression as the one reported by Cresser et al [s^ and, 
consequently, the same result. 

Jentschura [JT] pointed out that Cresser et aVs proce- 
dure is not gauge invariant since in a second order evalua- 
tion the sum over the complete spectrum of intermediate 
states is required to have equivalence between two dif- 
ferent gauges (more details are given in appendix of Ref. 

M)- 

Chluba and Sunyaev [43| developed another method to 
isolate the nonresonant contribution. In their method, 
the sum over all the intermediate states is split up in res- 
onant and nonresonant states. Although one can make 
conclusions for the difference between a pure cascade 
process, i.e., considering only the resonant states with 
a Lorentzian profile, and the two-photon emission given 
by all intermediate states (resonant and nonresonant), 
the definition of a nonresonant two-photon emission is 
unclear from a physical point of view. 

The values listed in Table IVII were used to calculate 
the nonresonant radiative corrections presented in Table 
1X1 (setting q = 10~'^). We notice that for the 3si/2 
1*1/2 transition the values calculated in this work differ 
from the values obtained by Jentschura [22j by 0.01 % 
for Z 1 and 0.1 % for Z = 40. 

The reason for some values in Table |X] being negative, 
such as the 3si/2 — *■ lsi/2 transition correction for Z = 
92, is due to the evaluation of the two-loop self-energy, 
which can be negative as any negative correction to the 
decay rates [22|. In Fig. [51 we represent the values of 
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the nonresonant radiative correction for several values of 
atomic number. 
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IV. CONCLUSIONS 



APPENDIX A: APPENDIX 



By applying a finite basis set constructed from B- 
splines to solve the Dirac equation, we have been able 
to calculate the decay rates in the Line Profile and QED 
based on the two-loop self energy approaches for all two- 
photon transitions from initial states with n = 2 and 3 
for a set of hydrogenlike ions with nuclear charge ranging 
from Z — 1 to Z = ^2. In these calculations the most 
significant multipoles contributions were considered, such 
as the 2EI, ElMl, 2MI, etc. We have also studied the 
spectral distributions of several transitions, which exhibit 
specific structures, such as resonances and transparen- 
cies. The latter reveal that two-photon emission is not 
possible at certain frequencies. The numerical results 
obtained in this work are in good agreement with other 
nonrelativistic and relativistic theoretical results. 

The QED approach gives a better contribution for 
a pure coherent nonresonant two-photon emission than 
Cresser's and Chluba's methods, not only because it is 
derived from physical arguments, but also due to the fact 
that it can be obtained from the Line Profile approach 
by removing the cascade process and setting the radia- 
tive corrections to zero. Therefore, it is a useful tech- 
nique in theoretical evaluations that require a coherent 
two-photon decay rate rather than the sum of this term 
along with the sequential one photon decay rate (cascade 
process). 

We conclude that the Line Profile approach is the most 
suitable for comparison with experimental results since it 
includes the terms associated with cascade process as well 
as radiative corrections. 

We end this conclusion by emphasizing that the 
method of integration used to obtain one electron de- 
cay rates (in both approaches and for both nonresonant 
and resonant transitions) can be adapted perfectly to ions 
with two or three electrons. 
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In order to make the task of deriving the matrix el- 
ements less cumbersome, further simplifications can be 
done in the matrix elements (Eq. (|14p and HI])) by notic- 
ing that the longitudinal part of the operator 'aL-^\,j [i^ , 



can be writen using a commutation relation as 



'la/ 



II 



D,(PL,M\ , 



(Al) 



(A2) 



where Ho stands for the Dirac Hamiltonian and (pL.M are 
the components of the spherical tensor of rank L resulting 
from the multipolc expansion of the potential A*. The 
reduction of Eq. (|A2|) to radial integrals along with the 
scalar term of the potential A*, lead to the following 

expression for the radial element matrix j ■* 



M 



i-l.L) 



G(2L + 1) 



-UJf^ 



(A3) 



where ujfi is the energy of the one- photon transition. 
This term is gauge independent for one-photon, as 
demonstrated by Grant [25j . since uifi — — w. Consid- 

2]), the radial matrix element m'^j'^^ can 



ering Eq. 

also be rewritten as 



M 



(2L+ 1) 



+ 

L-1 



LI 



L-l 



(A4) 



The explicit expressions of the derivatives of the matrix 
elements (Eq. Eq. and Eq. (|X3)) ') are given by 



d_ 



M 



2L + 1 




-^L{L + 2)J^' 



(A5) 



duj 

and 
d 

duj 



2L+1 



^^-1 " ^^^^ 
(A6) 



M 



Gi2L + l)U'^^^]4'-'' 



^[(i + l)c^ + (L + 2)c^/,] J(^)|, 

(A7) 
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where the integrals J^^^ and are defined by Using these expressions, along with the definitions 

given by Eqs. (|26p . we were able to obtain the coeffi- 

jR _ 1 / (P O + P O )rj (—^dr (A8) cients Oq and aj listed in Table IXll for the transition 
^ cJo ^ '^^'-^^{0) ' ^ ' 3si/2 lsi/2 for ions with Z = 1, 40 and 92. 



4"^^ ^ - r {PfP^ + QfQ^)r3L (—) dr. (A9) 
c Jo \ c / 
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TABLE I: Multipole contributions included in the present cal- 
culation of the total two-photon rate for the 2pi/2 ~^ ^^1/2 
transition. Comparison between the values obtained in this 
work (Eg. I17|l and other theoretical values. Powers of ten are 
given in parentheses. 



Multipoles 



Contribution (s ^) 





Z = 1 




Z = 40 


Z = 92 


ElMl 


9.676654(- 
9.667(-6) 
9.677(-6) 


-6) 

a 
b 


6.027323(7) 

6.020(7)" 

6.341(7)' 


3.863302(10) 

3.859(10)" 

4.966(10)' 


E1E2 


6.61179(- 
6.605(-6) 
6.673(-6) 


6) 

a 
h 


4.092020(7) 

4.088(7)" 

4.374(7)' 


2.358404(10) 
2.357(10) " 
3.425(10)' 


M1M2 


3.827877(- 


-17) 


5.602320(2) 


7.689142(6) 


E2M2 


9.385470(- 


-17) 


1.521687(3) 


2.834065(7) 


E2E3 


4.095985(- 


-18) 


6.608612(1) 


1.177403(6) 


Total 


1.628845(- 


-5) 


1.01195(8) 


6.225309(10) 



"Labzowsky et al 
''Labzowsky et al 



TABLE III: Transparencies for several two-photon transi- 
tions. The variable y = coi/iUfi is the fraction of the photon 
energy carried by one of the two-photons. 



Transition 

3si/2 lSl/2 

4Si/2 lSl/2 
6Si/2 lSl/2 



y{Z = l) y[Z = 40) y{Z = 92) 



0.780267 
0.7803 " 
0.7802 ' 
0.7803 " 
0.737322 
0.7373 " 
0.7373 " 
0.7032201 
0.7098 " 
0.7079 ' 
0.7098 " 



0.77628 



0.73273 



0.70497 



"Florescu et al [4J1 
''Quattropani et al [45 
■^Tung et al [sl 



0.7518 



0.7034 



0.6725 



TABLE II: Total two-photon decay rates (s~^) for the transi- 
tions 2si/2 —* lsi/2 and 2pi/2 — > lsi/2. Comparison between 
the values obtained in this work (Eq.[20} and other theoretical 
values. Powers of ten are given in parentheses. 









Total decay rate (s 




z = 


1 


/V 


2Sl/2 


2pi/2 






lSl/2 


8.229059 
8.2202 " 


1.628845(-5) 
1.6272(-5) " 
1.6350(-5)' 










z = 


40 


/V 


2si/2 


2pi/2 






lSl/2 


3.198851(10) 
3.1954(10) " 


1.01195(8) 

1.010(8)" 

1.071(8)' 


z = 


92 


./\* 


2Sl/2 


2pi/2 






lSl/2 


3.835978(12) 
3.8216(12) " 


6.225309(10) 

6.216(10)" 

8.391(10)' 



"Labzowsky et al 
''Labzowsky et al 



II 
3: 
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TABLE IV: Radiative corrections for several states in a.u. Powers of ten are given in parentheses. 



Radiative terms (a. u.) 



z = 


1 


state 


2Si/2 


2pi/2 


2p3/2 


3Si/2 


3pi/2 


3P3/2 






Re [SE]+VP= AS„ 
Im [SE]= rn/2 


1.5867(-7) 
1.9905(-16) 


-1.9542(-9) 
1.5162(-8) 


1.9095(-9) 
1.5162(-8) 


4.7376(-08) 
1.5281(-10) 


-5.5003(-10) 
4.5911(-9) 


6.4103(-10) 
4.5911(-9) 


z = 


40 


state 


2Si/2 


2pi/2 


2p3/2 


3Sl/2 


3pi/2 


3p3/2 






Re [SE]+VP= A£„ 
Im [SE]= r„/2 


8.6991(-2) 
1.4690(-6) 


-1.4913(-3) 
3.9208(-2) 


7.0905(-3) 
3.8037(-2) 


2.6401(-2) 
4.6373(-4) 


-1.3587-4) 
1.1689(-2) 


2.3058(-3) 
1.1627(-2) 


z = 


92 


state 


2Si/2 


2pi/2 


2p3/2 


3Si/2 


3pi/2 


3p3/2 






Re [SE]+VP= AS„ 
Im [SE]= r„/2 


1.7995 
4.7468(-3) 


2.4965(-l) 
1.1417 


3.2252(-l) 
9.5531(-1) 


5.7911(-01) 
2.7061(-2) 


9.3654(-2) 
3.1125(-1) 


1.102(-1) 
3.0444(-l) 



TABLE V: Sum of the terms \^^^ and f)i^*^, given by Eqs. (|28|l andl30|l. respectively, for transitions from bound states with 



Ui — 


3. 


This values were 


obtain using the radiative corrections of Table ITVl and q — 1. 


Powers of ten are j 


jiven in parentheses. 


Z = 


1 


/V 


3Sl/2 


3pi/2 


3p3/2 


3^3/2 


30(5/2 






lSl/2 


2.342758 


3.966941 (-6) 


3.230044(-6) 


3.706845 


3.706854 






2Sl/2 


6.452435(-2) 


4.925801(-8) 


4.926319(-8) 


7.762447(-4) 


7.750004(-4) 






2pi/2 


2.894796(-8) 


4.660148(-2) 


4.414498(-4) 


3.890718(-8) 


3.049476(-9) 






2p3/2 


5.789457(-8) 


8.832671 (-4) 


4.704893(-2) 


1.326766(-8) 


4.912187(-8) 


Z = 


40 


/V 


3Sl/2 


3Pl/2 


3p3/2 


3^3/2 


3^5/2 






lSl/2 


7.813052(9) 


2.872804(7) 


2.002016(7) 


1.436349(10) 


1.441446(10) 






2Sl/2 


2.247363(8) 


2.812583(5) 


3.395829(5) 


6.369235(6) 


-1.183514(6) 






2pi/2 


1.756667(5) 


1.797278(8) 


3.117290(6) 


3.012475(5) 


2.478360(4) 






2p3/2 


3.370367(5) 


8.545466(6) 


2.168435(8) 


8.878266(4) 


3.231816(5) 


Z = 


92 


/V 


3Sl/2 


3Pl/2 


3p3/2 


3d3/2 


3^5/2 






lSl/2 


-5.316586(10) 


3.647738(10) 


2.048248(11) 


1.646524(12) 


1.637676(12) 






2Si/2 


7.443297(9) 


1.683735(8) 


6.928218(8) 


1.135193(10) 


6.013487(9) 






2pi/2 


9.548939(7) 


2.314262(10) 


3.637390(9) 


4.841719(8) 


1.534768(7) 






2p3/2 


1.554405(8) 


5.312377(9) 


5.376493(10) 


1.660640(8) 


3.122434(8) 



TABLE VI: Sum of the terms i)'^^^ and i)J^^, given by Eqs. (|37[l and (|36|l . respectively, for transitions from bound states with 
Ui — 3. This values were obtain without radiative corrections, using q — 10~^ and following and using Jentschura's approach. 
Powers of ten are given in parentheses. 



z = 


1 


/V 


3Sl/2 


3pi/2 


3p3/2 


3^3/2 


30(5/2 






lsi/2 


2.342751 


3.966941(-6) 


3.230044(-6) 


3.706845 


3.706854 






2Sl/2 


6.452428(-2) 


4.925765(-8) 


4.926337(-8) 


7.762447(-4) 


7.750004(-4) 






2pi/2 


2.894793(-8) 


4.660148(-2) 


4.414498(-4) 


3.890718(-8) 


3.049476(-9) 






2p3/2 


5.789453(-8) 


8.832670(-4) 


4.704892(-2) 


1.326766(-8) 


4.912188(-8) 


z = 


40 


/V 


3Si/2 


3pi/2 


3p3/2 


3d3/2 


3^5/2 






lSl/2 


7.799875(9) 


2.872822(7) 


2.000453(7) 


1.436350(10) 


1.441394(10) 






2Sl/2 


2.246945(8) 


2.802618(5) 


3.400453(5) 


6.369529(6) 


-1.180529(6) 






2pi/2 


1.756241(5) 


1.797280(8) 


3.116223(6) 


3.080467(5) 


2.477876(4) 






2p3/2 


3.369973(5) 


8.514049(6) 


2.167575(8) 


8.872957(4) 


3.233897(5) 


z = 


92 


/V 


3Si/2 


3pi/2 


3p3/2 


3d3/2 


30(5/2 






lSl/2 


-5.843073(10) 


3.632557(10) 


2.037502 (11) 


1.648325(12) 


1.637187(12) 






2Si/2 


7.703355(9) 


1.624996(8) 


7.340967(8) 


1.135843(10) 


6.400342(9) 






2pi/2 


9.552389(7) 


2.313307(10) 


3.630875(9) 


4.857978(8) 


1.5282890(7) 






2p3/2 


1.556356(8) 


5.270175(9) 


5.351656(10) 


1.655470(8) 


3.134378(8) 
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TABLE VII: Same as TableUfor the transition 2p3/2 -> Isi/a 
Powers of ten are given in parentheses. 



Muhipoles 




Contribution (s ^~ 


) 




Z = 1 




Z = 40 


Z = 92 


ElMl 


9.700994(- 


-6) 


3.547078 (8) 


2.93820 (12) 


E1E2 


6.612242(- 


-6) 


4.597372(7) 


1.113120 (11) 


E1M3 


1.761532(- 


-16) 


3.264236(3) 


1.216566 (8) 


M1M2 


2.450145(- 


-15) 


4.265433(4) 


1.129082 (9) 


E2M2 


7.227055(- 


-17) 


1.337970(3) 


4.923824 (7) 


E2E3 


4.096369(- 


-18) 


7.718214(1) 


1.216566 (8) 


Total 


1.631323(- 


-5) 


4.007255 (8) 


3.050699(12) 



TABLE VIII: Same as Table |T] for the transition 3si/2 
2si/2. Powers of ten are given in parentheses. 



Multipoles 




Contribution (s ^) 






Z = 1 




Z = 40 


Z = 92 


2E1 


6.452436(- 


-2) 


3.167606(8) 


1.850546 (12) 


E1M2 


6.725935(- 


-14) 


1.150220(3) 


3.359725 (8) 


2M1 


1.038556(- 


-14) 


1.241685(2) 


1.031596 (6) 


2E2 


1.456030(- 


-14) 


1.584123(2) 


8.482772 (5) 


2M2 


1.901242(- 


-27) 


6.068107(^5) 


4.479357 (3) 



Total 6.452436(-2) 3.167620(8) 1.850884(12) 
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TABLE IX: Total two-photon decay rates (s ^) in the LPA, given by Eq. fEg. I20p . for transitions from bound states with 
Ui = 3. Powers of ten are given in parentheses. 



Total decay rate (s ^) 



z = 


1 


/V 


3Sl/2 


3pi/2 


3p3/2 


3^3/2 


3^5/2 








6.382020(6) 


3.431055(1) 


3.431043(1) 


7.213121(7) 


7.212970(7) 






2Sl/2 


6.452436(-2) 


4.925806(-8) 


4.965339(-8) 


7.762774(-4) 


7.751032(-4) 








6.4527(-2)'' 






7.7589(-4)" 








2pi/2 


2.894796(-8) 


4.6601485(-2) 


4.414514(-4) 


3.890719(-8) 


3.070354(-9) 










4.7484(-2)" 












2fl3/2 


5.789457(-8) 


8.832671(-4) 


4.704893(-2) 


1.326768(-8) 


4.912603(-8) 


z = 


40 


/V 


3Si/2 


3pi/2 


3p3/2 


3d3/2 


3^5/2 






lSl/2 


1.940069(13) 


1.465902(11) 


1.459393(13) 


1.909365(14) 


1.846466(14) 






2Sl/2 


3.167620(8) 


1.249212(6) 


6.204266(6) 


4.239616(8) 


1.256275(9) 






2pi/2 


2.488771(5) 


1.797286(8) 


2.342989(7) 


4.767668(5) 


1.268419(6) 






2p3/2 


3.370367(5) 


8.545466(6) 


2.648671(8) 


4.278490(5) 


3.734221(5) 


z = 


92 


/V 


3Sl/2 


3pi/2 


3p3/2 


3^3/2 


3^5/2 






lSl/2 


1.067206(15) 


5.649957(13) 


4.187334(13) 


6.170663(15) 


5.175621(15) 






2Si/2 


1.850884(12) 


4.581582(10) 


1.219838(11) 


5.062691(12) 


1.123469(13) 






2pi/2 


8.893987(9) 


2.336202(10) 


2.576780(11) 


9.697069(9) 


5.543826(10) 






2p3/2 


1.554405(8) 


5.312377(9) 


1.059702(12) 


1.959473(10) 


1.239531(9) 



"Tung et al iMl 



TABLE X: Total nonresonant two-photon correction (s^^) in the TLA , given by Eq. (|20|l . for transitions from bound states 
with Hi — 3. Comparison between the values obtained in this work and other theoretical values. Powers of ten are given in 
parentheses. 



Total nonresonant correction (s ^) 



z = 


1 


/V 


3Sl/2 


3pi/2 


3P3/2 


3d3/2 


3^5/2 






lSl/2 


2.082562 


2.981766(-6) 


2.98676(-6) 


1.042768 


1.042835 






2.082853 " 






1.042896 ' 








2Sl/2 


6.452428(-2) 


4.925721(-8) 


4.926293(-8) 


7.762407(-4) 


7.749962(-4) 






6.4530(-2) " 














2pi/2 


2.894793(-8) 


4.6601486(-2) 


4.414498(-4) 


3.890718(-8) 


3.049476(-9) 






2p3/2 


5.789453(-8) 


8.832670(-4) 


4.704892(-2) 


1.326767(-8) 


4.912188(-8) 


z = 


40 


/V 


3Si/2 


3Pl/2 


3p3/2 


3d3/2 


3^5/2 






lSl/2 


6.560351(9) 


2.224659(7) 


1.885681(7) 


3.456276(9) 


3.874677(9) 






lSl/2 


6.554(9) " 














2Sl/2 


2.245669(8) 


2.793230(5) 


3.395425(5) 


5.920457(6) 


-2.764917(6) 






2pi/2 


1.755215(5) 


1.797280(8) 


3.088316(6) 


3.078375(5) 


2.349043(4) 






2fl3/2 


3.369973(5) 


8.514049(6) 


2.166915(8) 


8.842254(4) 


3.233689(5) 


z = 


92 


/V 


3Si/2 


3pi/2 


3p3/2 


3d3/2 


3^5/2 






lSl/2 


-3.842113(11) 


2.891626(10) 


7.976296(10) 


8.916260(10) 


3.271857(11) 






2Sl/2 


5.570205(9) 


1.258571(8) 


7.132981(8) 


7.613467(9) 


-5.085413(9) 






2pi/2 


8.535648(7) 


2.313315(10) 


3.386706(9) 


4.767493(8) 


-5.866950(7) 






2p3/2 


1.556356(8) 


5.270175(9) 


5.253150(10) 


1.539227(8) 


3.127867(8) 



"Jentschura 25 
''Jentschura 46 
'^Jentschura 22 
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TABLE XI: Values of the coefficients (s"^) and (s^^), 
given by Eqs. (|26|l . of the transition 3si/2 —> lsi/2 for several 
values of Z. Powers of ten are given in parentheses. 





Z = 1 


Z = 40 


Z = 92 


1/2 


5.081524(-3) 


3.554138(10) 


3.793201(13) 


1/2 

a/ 


-3.468385(-l) 


-1.425932(9) 


-2.066515(11) 


3/2 

"o 


1.016398(-2) 


8.148532(10) 


1.250476(14) 


3/2 


-6.936930(-l) 


-2.929768(9) 


-3.4880108(11) 
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FIG. 1: Spectral distribution function il'{y, Z), defined by Eq. 
p9)) . of the ElMl and E1E2 contributions for the transition 
2pi/2 — > lsi/2 at Z=40. The variable y = Lo/uifi is the frac- 
tion of the photon energy carried by one of the two-photons. 
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FIG. 2: Spectral distribution function ip{y, Z), defined by Eq. 
(|39ll . for the transition 2pi/2 Isi/a at Z=l, 40, and 92. The 
variable y = u/uifi is the fraction of the photon energy carried 
by one of the two-photons. 




FIG. 3: Spectral distribution function tp{y, Z), defined by Eq. 
(|39|) . of the 2E1 contribution for the transition 3si/2 —> lsi/2 
at Z=l, 40, and 92. The variable y = u/uifi is the fraction of 
the photon energy carried by one of the two-photons. 
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FIG. 4: Transparency frequency, y^^^^^P^ of the transition 
3si/2 — * lsi/2 as function of the atomic number Z. 



1E10- 




>'=cd/co 



FIG. 5: Same as Fig. |3]for the transition 3d3/2 —> lsi/2- 
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_y=co/ft) 



FIG. 6: Spectral distribution function il>{y,Z), defined by 
Eq. (|39p . of the ElMl and contribution for the transition 
2p3/2 lsi/2 at Z=l, Z=40, Z=92. The variable y = uj/cofi 
is the fraction of energy carried by one of the two-photons. 



20 



1E-3 3 



1E-4 



1E-5 




^ 1 E-6 ^ 



1E-7- 



E\Ml 2p - 




ipTLA 

E\Ml 2p^,^- 




hr (^"') 






ElMl 2p -is 

^ Ml 1/2 





1E-8- 



— I — 
20 



— I — 
40 



60 



80 



FIG. 7: Multipole combination ElMl decay rate values 
Weimi, obtained in the LPA and TLA, as function of the 
atomic number Z for the transitions 2pi/2 lsi/2 (dash- 
dot) and 2p3/2 lsi/2 (solid and dot lines). The cascade 
term in LPA is represented by the dash line. 
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FIG. 8: Nonrosonant correction of the nmltipole combination 
2E1 in the 3si/2 — > lsi/2 transition divided by Z'^ as function 
of the atomic number Z. 



